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I present SUSY_FLAVOR - a Fortran 77 program able to calculate simultaneously 29 low-energy 
flavor and CP-violating observables in the general i?-parity conserving MSSM, including the 
case of large flavor violation in the sfermion sector. SUSY_FLAVOR v2 performs also the resum- 
mation of chirally enhanced corrections, arising in large tan/3 regime and/or large trilinear 
soft mixing terms, to all orders of perturbation theory. I discuss an example of application of 
SUSY_FLAVOR to analysis of the B° -> ^ + fj,~ decay rate in the MSSM. 



1 Introduction 



Flavor physics was in the recent years one of the most active fields in the high energy physics. 
New experiments helped to improve the accuracy of various measurements related to rare de- 
cays and put strong constraints on the flavor structure of physics beyond the Standard Model 
(SM), in particular imposing stringent limits on the flavor- and CP- violating parameters of the 
Minimal Supersymmetric Standard Model (MSSM). It is then increasingly important to have 
an universal computational tool which helps to compare new data with the predictions of the 
MSSM. Developing such a tool is a non-trivial task requiring extensive and tedious calculations. 
Numerous existing analyses usually consider only a few rare decays simultaneously and in most 
cases they are restricted to the case of so-called Minimal Flavor Violation (MFV) scenario, where 
the CKM matrix is the only source of CP and flavour violation^. 

Based on a series of papers where many rare processes were analyzed within the general 
MSSM a library of relevant computer codes has been published as SUSY_FLAVOR vl^. 
SUSY_FLAVOR v2^ is in addition capable of resumming leadi ng ch irally enhanced corrections, 
important in the regime of large tan f} or large trilinear Xterms™, to all orders of perturbation 
theory and for any pattern of sfermion mass matrices - a unique feature not shared by other 
publicly available programs. In this article I briefly summarize the main features of SUSY_FLAVOR 
and present example of its application to estimating the B® — > decay rate in the MSSM. 



2 SUSY_FLAVOR structure, input parameters and calculated observables 

SUSY_FLAVDR is capable of calculating physical observables within the most general i?-parity 
conserving MSSM, with one exception: currently it assumes massless neutrinos (and no right 
neutrino/sneutrino fields^), so the PMNS mixing matrix does not appear in the couplings. 

SUSY_FLAVOR has been in development long before the Les Houches Accord^] (SLHA) for 
common MSSM conventions was agreed on. Thus, the internal routines of the library follow 
the conventions of earlier paper^l. However, by default SUSY_FLAVOR can be initialized with 
a SLHA2 compatible set of parameters - all translations are done internally. Note that in 
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Table 1: One loop parton level formfactors implemented in SUSY_FLAVOR. /, J,K,L denote flavor indices. 



Box 


Penguin 


Self energy 


d l d J d K d L 
u I u J u K u L 
d z d J l K l L 
d I d J v K v L 


Zd L d J , jd'd J , gd L d-' 
H°d T d J , A\dfd J 
Hfu 1 ^, A9u T u J 
~fVl J 


d 1 -quark 
u / -quark 
charged lepton I 1 



SUSY_FLAVDR one can also use so-called "non-analytic j4-terms" of the form A'/" 7 H 2 * L 1 E J + 
A' d IJ H 2 *Q J D J + A'J J H^QtU 3 , which are not included in the standard SLHA2 parametrization. 
Calculations in SUSY_FLAVOR take the following steps^ 

1. Parameter initialization. Users can adjust the basic SM parameters and initialize all (or the 
chosen subset of) Higgs sector parameters and supersymmetric soft masses and couplings (which 
must be specified at the SUSY scale). 

2. Calculation of the physical masses and the mixing angles. In the next stage, the eigenvalues 
of the mass matrices of all MSSM particles and their mixing matrices at tree level are calculated 
numerically, without any approximations. 

3. Resummation of the chirally enhanced effects. In the regime of large tan/3 and/or large 
trilinear SUSY breaking terms, large chirally enhanced corrections to Yukawa couplings and 
CKM matrix elements arise^. They are resummed to all orders of perturbation theory using the 
formalism developed in^l. The level of resummation (no resummation, approximate analytical 
resummation in the decoupling limit, iterative numerical resummation) is a user defined option. 
4- Calculation of the Wilson coefficients at the SUSY scale. In the current version, SUSY_FLAVOR 
calculates Wilson coefficients generated by the diagrams listed in Table [1] (routines for given 
formfactor accept fermion generation indices as input parameters). 

5. Strong corrections. After evaluating virtual SUSY contributions, SUSY_FLAVOR performs the 
QCD evolution of the Wilson coefficients from the high (SUSY or top quark mass) scale to the 
low energy scale appropriate for a given decay. Necessary hadronic matrix element estimates 
and other QCD related quantities are treated as external parameters, initialized to the default 
values extracted from analyses done within the SM but also directly modifiable by users. 

6. Evaluation of physical observables. Finally physical observables are calculated and printed 
out. Current list of processes implemented in SUSY_FLAVOR v2 is listed in Table [2j 



3 Application of SUSY_FLAVOR to analysis of the B — > (j,*fx decay rate. 

One of the most promising signals for new physics at the LHC is the rare decay B® — > 
It is suppressed as a loop-level flavour-changing neutral current and by a lepton mass insertion 
required for the final state muon helicities. The LHC will be the first experiment able to probe 
this decay channel down to its SM-predicted branching ratio. The winter 2012 experimental 
95% CL bounds^ and the SM prediction^ for B® — > fJ. + fi~ decay rate can be summarized as: 

CMS < 7.7 x 10" 9 

LHCb < 4.5 x 1CT 9 

SM Prediction (3.35 ± 0.32) x 10" 9 
ATLAS, CMS and LHCb will be able soon to reconstruct the SM-like B® — > fi + [J>~ signal with 
significance of 3a, so that this very rare decay could be finally discovered and measured. 

SUSY_FLAVDR is an efficient tool allowing to estimate the size of possible SUSY contributions 
to the Bg — > n~ channel. In the MSSM, even in the restricted MFV case, for large values of 
tan/3 the B(B^ — > can be strongly enhanced^ {Ma is the CP-odd Higgs boson mass): 



Table 2: List of observables calculated by SUSY_FLAVOR v2 and their currently measured values. 



Observable 


Experiment 


Observable 


Experiment 


AF = 1 


AF = 




< 2 8 x If) -11 


?(S-2)e 
J(5-2)r 


(1159652188.4±4.3) x 10~ 12 


B(t — y e^) 


< 3 3 x 1CT 8 


(11659208.7 ±8.7) x lO" 10 


B(t — > u^y) 


< 4.4 x 10~ 8 


< 1.1 x 10~ 3 


B(Kt -± -K°VV\ 


< 6.7 x 10~ 8 


\d e (ecm) 


< 1.6 x 10~ 27 


B(K+ ->• 7T+W) 


!7.3±^;| x 10- 11 

< i 13 x irr 7 


d M |(ecm) 
|d T |(ecm) 


< 2.8 x 10~ 19 

< 1.1 x 10~ 17 




< 0.8 x 10~ 9 


d„|(ecm) 


< 2.9 x 10~ 26 


-)■ tt) 


< 4.1 x 10~ 3 


AF = 2 


B(B S ->■ ee) 
S(.B S -)■ /i/i) 


< 7.0 x 10~ 5 

< 4.2 x 10~ 9 


AM K 


(2.229 ±0.010) x 10" 3 
(5.292 ±0.009) x 10~ 3 ps" 1 


i3(B s -> tt) 




AM d 


(2.3718:??) x 10~ 2 ps- 1 


£(B S -> /ie) 


< 2.0 x 10~ 7 


AMs d 


(0.507 ±0.005) ps" 1 


i3(£ s -> re) 


< 2.8 x 10~ 5 


AM Bs 


(17.77 ±0.12) ps" 1 


i3(£ s -> pr) 

£(B+ -> T+^) 


< 2.2 x 10" 5 
(1.65 ±0.34) x 10~ 4 






B{B d ^Dlv) 

B(B -► X s7 ) 


(0.407 ±0.12 ±0.049) 
(3.52 ±0.25) x 10~ 4 







This result can be further significantly modified by non-vanishing flavor-violating terms in the 
sfermion mass matrices, leading to large contributions from box and Z-penguin diagrams. Apart 
from enhancing the decay rate, the interference of these terms could also conceivably lead to a 
cancellation that would suppress the branching ratio even below the SM prediction^. 

To quantitatively study the size of possible effects, one needs to perform a scan over the 
MSSM parameter space. The ranges of variation over MSSM parameters are shown in left panel 
of Table [3] (all parameters in scan are real; "SUSY scale" refers to the common mass parameter 
for the first two squark generations; the trilinear soft breaking terms are set to At = Af, = Mq^ 

and Af = M,-). Flavour violation is parametrized by the "mass insertions"^!, where I, J denote 
quark flavours, X, Y denote superfield chirality, and Q indicates the sfermion field: 



S&Y = 1 (M ^* y (2) 

l Q>YY 



(M*y x J x (M^ J - 



Realistic estimate of the allowed range for B(B® — > fi + n~) requires taking into account 
the experimental constraints from measurements of other rare decays. For that, all AF = 2 
observables, B(B — > X s j), B(Ki — y T[ a vv),B(K + — y -k + vu) decay rates and the neutron and 
electron electric dipole moments have been used out of quantities listed in Table [2J In addition, 
bounds on SUSY particle masses listed in right panel of Table [3] have been used. 

For the chosen bounds from Table [2] for which the experimental result and its error are 
known, parameter point in scan was accepted if 

\Qexp -Q th \< 3AQ exp + q\Q th \. (3) 

For the quantities for which only the upper bound is known, 

(1 + q)\Q th \ < Q exp (4) 



was required. The first and second terms on the right-hand side of Eq. [3] represent the 3a 
experimental error and the theoretical error respectively. The latter differs from quantity to 



Table 3: The range of input parameters and experimental constraints used for the numerical scan. All mass 

parameters are in GeV. 



Parameter 


Min 


Max 


Step 


tan/3 


2 


30 


varied 


CKM phase 7 





7T 


/or 

7T/25 


/"i 1 "1 1 1 T T * 71 f 

CP-odd Higgs Ma 


100 


500 


200 


Higgs mixing fi 


-40U 




oUU 


K\1C)\ winn mass A/f n 


1UU 


^00 




Gluino mass M3 


3M 2 


3M 2 





SUSY scale M SUSY 


500 


1000 


500 


Slepton Masses Mjr 


Msusy 


Msusy 





Left stop Mq l 


200 


500 


300 


Right sbottom Afr 

Or 


200 


500 


300 


Right stop Mi 


150 


300 


150 


dLL'' "dLL 


-1 


1 


0.1 


£13 X23 


-0.1 


0.1 


0.01 



Mass 



m { 
mh 



Constraint 



> 46 

> 94 

> 89 
> 95.7 



> 92.8-114 depending on sin 2 (a — f3) 



'0 which was assume generically to account 
6 for details). 



quantity and is usually smaller than the value q 
for the limited density of a numerical scan (see ref.' 

Fig. [T] shows the predictions for B{B® — > h + ijl~) over a general scan of 20 million points 
according to Table 5^ LL (left panel) and S^ 3 LR (right panel) were varied one at a time while 
setting the other to zero, e.g. all 5^ Y = and only 5^ 3 LL 7^ in the left panel. When 5f , LL 
is varied in the range [—1,1], one finds B(B® — > n + fi~) m in ~ 10~ 9 . This minimum is almost 
independent of tan /3. |^^/J can take on values up to « 0.9 and still pass all imposed constraints, 
though points beyond 0.3 are less dense. More interesting is the case when S^ LR is varied in 
the range [—0.1,0.1]. One can find a narrow cancellation region around S^ 3 LR ~ —0.01 and 
tan/3 < 10 where B(B® — > fJ, + fJ,~) m i n ~ 10~ 12 . This is three orders of magnitude lower than the 
SM prediction, making it effectively unobservable at the LHC. 
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Figure 1: Predictions for B(B S — > ) vs tan/3 from the scan of MSSM parameters in Table [3] Left(right) 

panel: 8^ LL varied. The dashed line shows the SM expectation. 



4 Conclusions 

I have presented SUSY_FLAV0R, a tool capable of calculating simultaneously 29 important flavor 
observables in the general -R-parity conserving MSSM. The calculation of the SUSY tree-level 
particle spectrum and flavor mixing matrices are performed exactly. SUSY_FLAV0R v2 implements 
also the resummation of chirally enhanced corrections, stemming from large values of tan (3 
and/or large trilinear ^4-terms. Thus SUSY_FLAV0R v2 is valid for the whole parameter space 



of the general i?-parity conserving MSSM, without restrictions on the size of the off-diagonal 
elements in the sfermion mass matrices - a unique feature currently not shared by other publicly 
available programs calculating FCNC and CP violation in supersymmetric models. I hope that 
SUSY .JFLAVOR becomes an important tool useful both for theorists and experimentalists who need 
to perform multi-process flavor analyses within the MSSM. 

As an example of such analysis, SUSY_FLAVOR has been used to perform a numerical explo- 
ration of the MSSM parameter space and estimate of the B® — > /i + /U~ decay rate. Scan shows 
that there exist cancellation regions where the contribution of diagrams with supersymmetric 
particles interferes destructively with SM diagrams, thus allowing the branching ratio to be 
significantly smaller than the SM prediction. Such effects may effectively hide the dimuon B® 
decay mode from the LHCb even though it is supposed to be one of the experiment's benchmark 
modes. Barring such cancellations, supersymmetric contributions typically tend to enhance the 
branching ratio for B® — > fi + jjL~ even for moderate values of tan (5 < 10 so that an experimen- 
tal measurement close to the SM prediction puts strong bounds on the size of allowed flavour 
violation in the squark sector. 



SUSY_FLAVDR can be downloaded from the address http : / / www . f uw . edu . pi/ susy jflavor 
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